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Hybrid systems consisting of different types of qubits are promising for building quantum com¬ 
puters if they combine useful properties of their constituent qubits. However, they also pose addi¬ 
tional challenges if one type of qubits is more susceptible to environmental noise than the others. 

Dynamical decoupling can help to protect such systems by reducing the decoherence due to the 
environmental noise, but the protection must be designed such that it does not interfere with the 
control fields driving the logical operations. Here, we test such a protection scheme on a quantum 
register consisting of the electronic and nuclear spins of a nitrogen-vacancy center in diamond. The 
results show that processing is compatible with protection; The dephasing time was extended almost 
to the limit given by the longitudinal relaxation time of the electron spin. 
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-Introduction. Choosing the most suitable qubit systems 
and implementing high-precision quantum gate opera¬ 
tions on them may be considered as the main challenges 
for building quantum computers that will be capable of 
solving certain problems qualitatively faster than classi¬ 
cal computers 0, Q- In every physical realization, en¬ 
vironmental noise is present, which degrades the infor¬ 
mation stored in the quantum systems. Several tech¬ 
niques have been developed for protecting quantum sys¬ 
tems against this effect, including passive schemes like 
decoherence-free subspaces [Stl or active techniques like 
dynamical decoupling (DD) [J-Q- So far, these tech¬ 
niques have been applied mostly to the protection of 
quantum states in quantum memories |7Hl2j|. However, 
similar protection is also required for quantum informa¬ 
tion processing. Here, the protection must be designed in 
such a way that it decouples the protected system from 
environmental noise, but not from the control fields that 
are applied for driving the gate operations. Experimental 
implementations of protected single-qubit gates and some 
special 2-qubit gates, such as a controlled-NOT (CNOT) 
gate, were recently reported [l3l - [l^ . but so far no scheme 
was demonstrated that can be applied to arbitrary two- 
qubit gates in different physical qubit systems. 
Implementing quantum computing in hybrid systems is 
a promising strategy for building quantum computers, as 
it increases the range of possible systems and allows one 
to combine useful properties of different types of qubits 
(TtI - I^ . As a typical hybrid system, the nitrogen va¬ 
cancy (NV) center in diamond consists of an electronic 
spin 5" = 1 and a nitrogen nuclear spin 1=1 and com¬ 
bines several useful properties for a quantum register. It 
also poses some challenges, mostly because the character¬ 
istic properties of the two types of spins differ by many 
orders of magnitude. For example, the natural coherence 
time at room temperature of the electron spin is of the 
order of 1 ^s, while that of the nuclear spin is 5 ms [l5l| . 
The electron spin dephasing time is usually shorter than 
typical durations of gates applied to the nitrogen spin. 
The critical gate operations in such a system are there¬ 


fore the two-qubit gates in which the electron spin is the 
control quit, while the nuclear spin is the target qubit: In 
this situation, the relatively slow Rabi frequency of the 
nuclear spin results in a long duration of the gate, dur¬ 
ing which the whole quantum information must be con¬ 
served, including superpositions of different basis states 
of the electron spin. 

It is possible to apply DD pulses that act parallel to the 
control field driving the targeted nuclear spin [l^ . How¬ 
ever, since the DD pulses interchange the computational 
basis states of the control qubit (the electron spin), the 
segments of the control field separated by the DD pulses 
implement two gates with different control conditions. 
Obtaining the targeted evolution requires therefore care¬ 
ful synchronization of the DD pulses with the internal 
system evolution. In some systems, specific properties 
of the Hamiltonian allow one to design DD sequences 
such that they protect the system from environmental 
noise and simultaneously implement a useful gate oper¬ 
ation [l^. However, this approach is very specific to 
the system and requires that the quantum system can be 
simulated classically. It is therefore difficult to extend 
this approach to large quantum registers. 

Here, we present a strategy that avoids these limitations 
and can be applied to arbitrary gate operations with min¬ 
imal overhead. We demonstrate the scheme for the case 
of a 2-qubit system, where the duration of one gate is 
long enough that decoherence destroys a significant part 
of the quantum information, unless DD is applied to pro¬ 
tect it. Compared to existing schemes, our approach can 
implement arbitrary 2-qubit gates and does not require 
fine-tuning to system parameters. It requires that the 
control field driving the gate operation is split into serval 
segments that can be inserted into the delays of the DD 
sequence. The segments have to be adjusted to take the 
effect of the DD pulses into account, which interchange 
the computational basis states |0) e->- |1) and therefore 
the condition for the controlled 2-qubit gate. 
-Experimental protocol and results. We illustrate this 
basic idea for the case of a controlled rotation, CR. The 
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FIG. 1: Circuits for the controlled-rotation gate and its pro¬ 
tected version. Rx{0) represents a rotation around the x-axis 
by an angle 9. The rotation angles are related by 261+62 = 9. 
The filled rectangles are the DD pulses with flip angle tt. 


unitary operation for this gate can be represented as 

= (I i!°(e))' 

where 1, 0 and and Rxifi) represent 2x2 matrices corre¬ 
sponding to the unit operator, the zero matrix and a ro¬ 
tation matrix around the x-axis, Rx{0) = where 

ax denotes the x-component of the Pauli matrix. 

FigH] shows the circuit models of the CR gates. The 
left-hand part represents the unprotected version, while 
the right-hand side represents the protected gate for the 
simplest DD sequence consisting of two refocusing pulses. 
Filled circles indicate that the operation is performed on 
the subspace where the control qubit is in the logical 
state |1), empty circles stand for the control condition 
|0). When 9 = tt, Ucr becomes the CNOT gate, with an 
additional phase gate on qubit 1 and an overall 

phase. 

For the experimental implementation, we chose a 4-level 
system out of the 9-level system composed by the elec¬ 
tronic and nuclear spins of the NV center in diamond 
shown in Fig. [21 If the magnetic field is oriented along 
the NV symmetry axis, the Hamiltonian of the NV center 
can be written as [2lj 

Hs = DSl + + Pll + -fnBh + ASJ, ■ ( 2 ) 

Here Sz and Iz are the 2 -components of the spin-1 oper¬ 
ators for the electron and nitrogen (^'^N) nuclear spins, 
respectively. The zero-field splitting is D = 2.87 GHz, 
the nuclear quadrupolar splitting is P = —4.95 MHz, 
and the hyperhne coupling is H = 2.16 MHz [H, [^ . 
The electronic gyromagnetic ratio is 7 e = 2.8 MHz/G, 
and the nuclear gyromagnetic ratio 7 „ = 0.30 kHz/G. 
The strength of the field in our experiment was about 87 
G. 

In the 9-dimensional Hilbert space of the electronic and 
nuclear spins, we chose two different 2-qubit systems to 
implement the gates. In Fig. |21 they are denoted as a 
and 6, respectively. In each system, we use the electron 


DS^^ + YeBSz + YnBI z + AS^D 



FIG. 2: (color online). Structure and energy level scheme of 
the NV center. Two subspaces of the total system, which are 
used in the experiment, are shown and marked as subsystems 
a and b. The quantum numbers refer to the electronic spin S 
and the nuclear spin 7, respectively. 



Transition 
frequency (MHz) 

Rabi 

frequency (kHz) 

System a 

ui = 4.970 

U 2 = 2.808 

Di/( 27 r) = 9.05 
D 2 /( 27 r) = 5.84 

System b 

- 4.918 

U 2 = 7.088 

Di/( 27 r) = 9.00 
D 2 /( 27 r) = 5.60 


TABLE I: Measured parameters of the nuclear-spin transi¬ 
tions. 


spin as the control qubit (qubit 1), while the nuclear spin 
of the nitrogen represents the target qubit (qubit 2). The 
logical states for the two qubits correspond to 

{|0)i,|l)i} = {|0)e,|-l)e} 

{10)2,11)2} = {| 0 )„,|Tl)n}, ( 3 ) 

where the minus sign refers to system a and the plus sign 
to system b. Table U lists the measured transition fre¬ 
quencies of the nuclear spin and the corresponding Rabi 
frequencies. 

In the experiment, we used a enriched diamond as 
the sample, so that decoherence due to ^^G nuclear spins 
is small. The experiment was performed at room tem¬ 
perature. Fig. O shows the experimental procedure and 
the pulse sequence. Since the electron spin decoheres on 
a much shorter timescale than the nuclear spin, it is suffi¬ 
cient to apply DD to the electron spin transitions (qubit 
1). This is done in parallel to the controlled rotation 
of qubit 2. The controlled rotation of qubit 2 is inter¬ 
rupted during the application of the DD pulses, which 
are several orders of magnitude shorter than the dura¬ 
tion of the nuclear spin gate operations (10 ns vs. > 100 
/is). The DD pulses are hard 7r-pulses around the x- or 
y- axis, as indicated above the rectangles in Fig. [31 with 
a Rabi frequency of about 11 MHz. For this demonstra¬ 
tion experiment, we used only two DD pulses. To make 
it symmetric with respect to the initial condition, we al¬ 
ternated the rotation axes between x and y. To turn this 
into a cyclic operation (generating a unit operation), it 
has to be combined with a rotation around the 2 -axis, 
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FIG. 3: (color online). Pulse sequence and experimental pro¬ 
cedure. The microwave (MW) pulses are resonant with the 
transition |0)e|0)n -n- \ — l)e|0)„. In step 1) and 3), the MW 
pulses have the same Rabi frequency of about 0.28 MHz and 
the same flip angle. These pulses work as transition selective 
pulses. The DD pulses in step 2) are hard pulses with Rabi 
frequency of about 11 MHz. The frequencies of the radio 
frequency (RF) pulses are indicated in the rectangles, where 
v = U 2 , v' = I'l for a CR gate with control state |1) or = zzi, 
v' — U 2 for the gate with control state |0). 


which we implemented as a phase shift of the subsequent 
pulses. 

In step 1), we prepared the input state for the CR gate. 
The laser pulse initializes the electronic spin in its ground 
state nis = 0, while the nuclear spin is unpolarized. The 
first microwave (MW) pulse drives selectively the |00) -O- 
|10) transition, with a 7r/2 flip angle around the x-axis. 
The resulting input state is 

P^n = i(|00) - *|10))((00| + z(10|) + i|01)(01|, (4) 

written in the computation basis ([3]) but omitting the 
index. In this state, the control qubit is in a superposition 
of |0) and |1). The CR gate drives the transition labelled 
vi in Fig. [^if the control condition is 0 and the transition 
labelled V 2 if the control condition is 1. 

For the readout, we measure the population of the |0) 
state of the electron spin. To measure the effect of the 
RF pulse, which affects only the nuclear spin, we there¬ 
fore need to transfer the information from the nuclear to 
the electron spin system. For this purpose, we use again 
a selective MW pulse that rotates the |00) O |10) transi¬ 
tion around the x-axis, by a flip angle —a. As a function 
of the rotation angle 6 of the CR gate, the expected sig¬ 
nal (i.e. the occupation probability for the |0)-state in 
the computational space) is 

s= ^[l + cos^]^ + ^ = g[l + 2cos^-bcos^^]-bi (5) 

for an ideal gate. In the case of an unprotected gate, 
the environmental perturbation causes dephasing of the 
electronic spin state. For a Markovian environment, the 
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FIG. 4: (color online). Experiment results for implementation 
of the CR gates starting with the input state pin. The results 
for the unprotected gate in system a and the protected gate 
in systems a and b are shown in figure (a-c), respectively. The 
theoretical expectations for ideal gates are shown as red dash- 
dotted curves, and the experimental data as empty circles. 
The filled circles show the decay of the input state, without 
CR gate and without DD pulses in figure (a), and with the 
DD pulses in hgures (b-c). The dashed curves show the result 
of a linear fit. In hgure (a), the solid curve shows the result 
of a simulation with the measured dephasing time. In figure 
(b-c) the solid curves show the result of a fit to the function 
Sj. = (1 — Kt)sT 2 , where k is an empirical decay rate. 


off-diagonal elements of the electron spin qubit decay as 
g-t/'T 2 ^ This changes the detected signal to 

Ir & -f/T 2 0, 1 

st 2 = g[l+ 2cos-e -bcos -]-t--. (6) 

We implemented the unprotected gate by setting v' = v 
and switching off the DD pulses in the sequence in Fig. [S] 
Fig. 111(a) illustrates the resulting signal when the unpro¬ 
tected CR gate was implemented with control state |0) 
in system a. The horizontal axis gives the rotation angle 
of the CR gate 9 = flit, where Di is the Rabi frequency 
of the nuclear spin transition. While the ideal gate gen¬ 
erates a signal that is periodic in 6 with a period of 47r, 
the dephasing of the electron spin eliminates this signal 
component over a time scale of the order of the decay 
time T 2 ~ 34 ^s of the input state without DD pulses. 
The remaining signal component is periodic with 27r, in 
agreement with Eq. ([6]). The measured data for the CR 
gate are marked as empty circles, and the theoretical ex¬ 
pectation s for the ideal gate is shown as the dash-dotted 
curve. The solid curve shows the function st^ from Eq. 
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(|6|) and is in good agreement with the experimental data. 
Clearly, the dephasing time in this system is too short to 
allow the implementation of a nontrivial gate, such as 
CNOT, with sufficient precision. In addition, we mea¬ 
sured the decay of the initial state pin due to relaxation, 
without applying a gate operation. The resulting data 
points are shown as filled green circles. 

Figs, m (b-c), show the corresponding data obtained by 
the protected CR gate in systems a and b. The total 
rotation angle is 0 = (fli -I- n2)t/2. To evaluate the per¬ 
formance of the gate, we first measured the decay of the 
input state pin protected by the DD pulses. The results 
are shown in Fig. |4] (b-c) as filled green circles. Over the 
measured time scale, the experimental data points decay 
with a time constant of 1 /k ss 2.4 ms, mainly from lon¬ 
gitudinal relaxation (Ti) and experimental errors. For 
the protected gate, we used the same decay function and 
superimposed it over the experimental data points in fig¬ 
ures (b-c). The experimental data points can then be 
fitted with Sr = (1 — Kt)sT 2 , even for times much longer 
than the dephasing time scale in the unprotected gate 
shown in figure Figs. 0] (a), using pure depahsing times 
T 2 of 4 ms. 

A precise check of the performance of the gate operation 
could be performed by quantum process tomography. 
However, this requires the application of SWAP oper¬ 
ations consisting of the CR gates that we are evaluating. 
We therefore only checked the state of the control qubit 
after the gate operation by quantum state tomography. 
The CR gate is an entangling gate. However, when the 
rotation angle 9 is an even multiple of 27r, the CR gate 
0 becomes the identity or phase gate for qubit 

1, the electronic qubit can be separated from the nuclear 
spin. At these points, we can reconstruct the state of 
the electronic qubit by treating it as a single qubit. We 
implemented quantum state tomography for 0 = 0, 27r 
and 47r. To reconstruct each state, we used five measure¬ 
ments with different MW pulses in step 3): no pulse, and 
pulse along x, —x, y and —y, respectively. Fig. [S] illus¬ 
trates the results. One can clearly observe the effect of 
the 2 -rotation i.e., phase flip, in the results for 9 = 2tt. 
We compared the experimentally measured states with 
the ideal states by calculating the fidelity F = 
\Tr{pexpPideai}\ and obtained the values F = 0.90, 0.66 
and 0.57 for 0 = 0, 27r and 47r, respectively. We can 
attribute the difference between the experimental and 


ideal results mostly to longitudinal relaxation, which was 
measured to be « 3.5 ms by the standard MW pulse se¬ 
quences and cannot be reversed by dynamical decoupling, 
and to experimental imperfections. 

In addition to the experiments discussed here, we applied 
the protected CR gate also to other input states and 
found good agreement with the theoretical predictions. 
In those cases, where the initial state does not contain 
coherence of the electron ^n, the unprotected gate also 
works relatively well . The resulting evolution is 

described by Eq. m with T 2 —?> 0, and it corresponds to 
the results shown in Fig. 0] (a) for 9 > 2tt. 

0 = 0 0 = 2n 0 = 4n 



FIG. 5: (color online). Results for state tomography for 0 = 0, 
2n and drr shown as columns 1-3, respectively, where the first 
and second rows show the real and imaginary parts. In theory, 
density matrices are {E-\-ay)/2, [E — Gy) (2 and {E + Gy) (2 for 
0 = 0, 27r and dyr, respectively, where E denotes the identity 
operator. 


-Conclusion. We have demonstrated the protection of 2- 
qubit gate operations by dynamical decoupling by com¬ 
bining the control field with dynamical decoupling of the 
control qubit, which otherwise undergoes rapid dephas¬ 
ing. The protocol can be generalized for multiple-qubit 
gates, e.g., to NV centers in natural abundance diamond 
samples, where nuclear spins can act as qubits. The 
same strategy and techniques should be helpful for im¬ 
plementing quantum control in other quantum systems, 
in particular for hybrid systems. In future work, we plan 
to generalize this scheme and use more advanced DD se¬ 
quences, such as the KDD sequence (25j |. 
Acknowledgments.- JZ acknowledges F. D. Brandao for 
helpful discussions. This work is supported the DFG 
through Su 192/19-2. 


[1] M. Nielsen and I. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[2] J. Stolze and D. Suter, Quantum Computing: A Short 
Course from Theory to Experiment (Wiley-VCH, Berlin, 
2nd edition, 2008). 

[31 D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 
Lett. 81, 2594 (1998). 

[4] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 


[5] A. M. Souza, G. A. Alvarez, and D. Suter, Phil. Trans. 
R. Soc. A 370, 4748 (2012). 

[6] P. Zanardi, Physics Letters A 258, 77 (1999) 

[7] J. Du, X. Rong, N. Zhao, Y. Wang, J. Yang, and R. B. 
Liu, Nature 421, 1265 (2009). 

[8] G. deLange, Z. H. Wang, D. Riste, V. V. Dobrovitski, 
and R. Hanson, Science 330, 60 (2010). 

[9] G. A. Ryan, J. S. Hodges, and D. G. Gory, Phys. Rev. 
Lett. 105, 200402 (2010). 







5 


[10] Y. Wang, X. Rong, R Feng, W. Xu, B. Chong, J.-H. Su, 
J. Gong, and J. Du, Phys. Rev. Lett. 106, 040501 (2011). 

[11] C. Zu, W.-B. Wang, L. He, W.-G. Zhang, C.-Y. Dai, F. 
Wang, and L.-M. Duan, Nature 514, 72 (2014). 

[12] J. Shim, I. Niemeyer, J. Zhang, and D. Suter, Eur. Phys. 
Lett. 99, 40004 (2012). 

[13] A. M. Souza, G. A. Alvarez, and D. Suter, Phys. Rev. A 
86, 050301 (2012). 

[14] J. Zhang, A. M. Souza, F. D. Brandao, and D. Suter, 
Phys. Rev. Lett. 112, 050502 (2014). 

[15] T. van der Sar, Z. H.Wang, M. S. Blok, H. Bernien, T. H. 
Taminiau, D. M. Toyli, D. A. Lidar, D. D. Awschalom, R. 
Hanson, and V. V. Dobrovitski, Nature 484, 82 (2012). 

[16] G.-Q. Liu, H. G. Po, J. Du, R.-B. Liu, and X.-Y. Pan, 
Nat. Commun. 4, 2254 (2013). 

[17] H. Wu, R. E. George, J. H. Wesenberg, K. M0lmer, D. 1. 
Schuster, R. J. Schoelkopf, K. M. Itoh, A. Arda- van, J. 
J. L. Morton, and G. A. D. Briggs, Phys. Rev. Lett. 105, 
140503 (2010). 

[18] Y. Kubo, F. R. Ong, P. Bertet, D. Vion, V. Jacques, D. 
Zheng, A. Dreau, J.-F. Roch, A. Auffeves, F. Jelezko, et 


ah, Phys. Rev. Lett. 105, 140502 (2010). 

[19] D. 1. Schuster, A. P. Sears, E. Ginossar, L. DiCarlo, L. 
Frunzio, J. J. L. Morton, H. Wu, G. A. D. Briggs, B. B. 
Buckley, D. D. Awschalom, et ah, Phys. Rev. Lett. 105, 
140501 (2010). 

[20] 1. Ghiorescu, N. Groll, S. Bertaina, T. Mori, and S. 
Miyashita, Phys. Rev. B 82, 024413 (2010). 

[21] Y. Wang, F. Dolde, J. Biamonte, R. Babbush, V. 
Bergholm, S. Yang, 1. Jakobi, P. Neumann, A. Aspuru- 
Guzik, J. D.Whitfield, et ah, arXiv: 1405.2696 (2014). 

[22] G. S. Shin, M. G. Butler, H.-J. Wang, C. E. Avalos, S. J. 
Seltzer, R.-B. Liu, A. Pines, and V. S. Bajaj, Phys. Rev. 
B 89, 205202 (2014). 

[23] X.-F. He, N. B. Manson, and P. T. H. Fisk, Phys. Rev. 
B 47, 8816 (1993). 

[24] J. H. Shim, 1. Niemeyer, J. Zhang, and D. Suter, Phys. 
Rev. A 87, 012301 (2013). 

[25] A. M. Souza, G. A. Alvarez, and D. Suter, Phys. Rev. 
Lett. 106, 240501 (2011). 


